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Abstract  of  Thesis  Presented  to  the  Graduate  Council 

in  Partial  Fulfillment  of  the  Requirements  for  the 

Degree  of  Master  of  Science  in  Engineering 

THE  EFFECT  OF  CURVATURE  ON  THE  TRANSITION 
TO  TURBULENCE  IN  CIRCULAR  PIPES 

By 
Addison  Guy  Hardee,  Jr. 
December,  19&5 

Chairman:   Dr.  K.  T.  Mlllsaps 

Major  Department:   Aerospace  Engineering 

The  experimental  determination  of  the  critical 
Reynolds  number  for  transition  from  laminar  to  turbulent 
floxv  in  a  curved  pipe  of  circular  cross  section  has  been 
accomplished  by  the  use  of  a  streaming  potential  tech- 
nique.  The  streaming  potential  was  generated  by  an 
electrolyte  (0.001  H  KC1)  flowing  through  a  capillary 
tube  with  a  diameter  of  0.0709  inches  with  a  radius  of 
curvature  of  55«9  inches.   At  the  onset  of  turbulence, 
the  streaming  potential  began  to  fluctuate  in  an  early 
identifiable  manner.   This  fluctuation,  combined  with 
mass  flow  measurements,  was  used  to  determine  the  critical 
Reynolds  number  to  be  3190  1  10  for  a  fixed  geometry  and 
fluid. 
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I.   INTRODUCTION 

1.   Scope 

> 

The  purpose  of  the  present  investigation  is  to 

determine 

the  effect  of  curvature  on  the  transition  from 

laminar  to 

turbulent  flow  in  pipes  of  circular  cross  sec 

tlon.   The 

technique  used  in  observing  the  onset  of  tur- 

bulence  is 

one  developed  by  Anderson^)  Involving  the 

electric  double  layer. 

2-  ^^^mesusmom   Concerning  TTOf-^tn  ,H 

1               The  experimental  Investigations  of  flow  In 

ourved  pipes  may  be  divided  Into  two  main  olasses:   those 

prior  to  a 

Paper  by  Dean(2),  and  those  followlng>   ^ 

showed  from 

a  mathematical  analysis  of  the  equations  of 

motion  that 

dynamical  similarity  depends  only  on  the 

parameter 

K'  =(&VWUX) 

D   \*)\     -v-  -J                                                       (1) 

where     a 

is  the  radius  of  the  pipe  bore. 

R 

is  the  radius  of  curvature  of  the  pipe. 

Wo 

is  the  mean  velocity  of  the  flow  through  a 

1 

straight  pipe  under  the  same  pressure 
gradient, 

and   v    is  the  kinematic  viscosity  of  the  fluid. 

Of  course,  dimensional  analysis  shows  that  any  dynamical 
quantity  depends  on  the  functional  form  £,("r  ,  ^r-)  » 
where  -*~-   is  a  Reynolds  number.   A  particularly  simple 
form  is  ^- 

To  do  this  Dean  considered  the  flow  of  an 
incompressible  fluid  in  the  coordinate  system  where  (r.H'.e) 
are  the  orthogonal  coordinates  (specified  in  Fig.  1).   The 
corresponding  velocity  components  are  (U,V,W)  respectively. 
It  is  assumed  that  the  flow  is  steady  and  that    U ,V 
and  W   ,  but  not  the  pressure,  are  independent  of  ©  . 

The  general  orthogonal  curvilinear  form  of  the 
Navier-Stokes  equations  and  the  continuity  equation  are 
given  for  reference  in  Appendix  I.   In  the  transformation 
to  Dean's  coordinate  system 

1  =■  R  -+  r  sin  ^  (2) 


"x=  "IcosS  =  (R-*rsin  4»)  cos  Q 

y  =  I  sin  e  ■  (R  +  rs:nV)  sin  0  (3) 

7.    -   r  Co?  & 


Fig.  1.   Dean's  coordinate  system  for  flow  in 
a  curved  pipe. 


< 


and  therefore  the  scale  factors  are 


k  -  WW^ 


^--mNMf+m 


W  h*\x     pi'* 


k=vu-w  n^, 


Thus  to  write  the  Navier-Stokes  equations  of  motion  for 
Dean's  coordinate  system  we  use  the  relations 


X,-* 

,     w,»  u 

xz  -4' 

,       **«V 

>}r    G 

,      u^W 

(5) 


K,m  1 


^t"  r  (6) 


H3  r  R  i-  rs  in  Q 


The  continuity  equation  becomes 


i^  t  ll      tJs-w  y    +.  -L  iy  +  ^c°^  V   -  o      (7) 


Denoting  the  vector  sum  of  the  velocity  com- 
ponents as  v"  and  since  d.v  v  -O,  we  may  add  the  ex- 
pression Bs^U'^v)   (where  B  is  any  constant)  to  the 
Navier-Stokes  equations.   To  the  first  equation  of 
motion  we  add  2V5^-(<J;v  vj  ,  to  the  second,  add  ^F  5y  U,v  W. 
Nothing  is  added  to  the  third.   The  equations  of  motion 
are  then 


iiiy.  +  Y  iU  _  Yl  _  wVjlI  _  _WP\ 
u  3<*    r    zy        t         R-^sinY      ^  V5"' 


V"  5f   R  +  nmvA^r   r      r    ly 

ii^+y^v     uv  _  wxcos »    _  _j_  s_/ip_\ 

u  ^   Mt  T  r    R  +  rsmt       r  » ^  [  f>  J 


*   \^<    R+  rs 


&V 


V 


R  +  rs«n  vy/l  ^r  +  r 


and 


°T  r     W  R+rsm^  R+rsih4> 


■+ 


(9) 


i  a  /  i  aw  +    Wco,  y    )\  do) 


With  the  assumption  that  Vr  is  small,  the 
equations  of  motion  become 


=  l_/£\     2i^^  +  y-l-^)  (ID 


i  j^y    y.  <^  -+.  yy  _  w  cos  v 


and 


-+   "V 


y\/V  ar  \^r         «"        <C    *4>, 


(12) 


and  the  continuity  equation  is 


8 
Since   U  ,  V   and   W  are  assumed  to  be  indepen- 
dent of  6  ,  then  from  (13)  it  may  be  seen  that  V/p>    must 
have  the  form 


eUr.HO  +  U^Y) 


From  (11)  and  (12)  it  may  also  be  seen  that  f^_  must  be  a 
constant.   Thus  we  can  write 


where   G  is  a  constant  and  is  the  mean  pressure  gradient. 
We  may  satisfy  the  continuity  equation  with 


-ii-     af         v-  ^  (16) 


where  f«  ^C*-,^)  .   Using  relations  (16)  in  the  equations  of 
motion  and  eliminating   p   from  (11)  and  (12),  we  have  the 
two  equations 

-  _  T/f^  f  (17) 


and 


(18) 


where 


\7*  =     -^-    +1L    +  i-  ^-  (19) 


We  may  put  equations    (17)    and   (18)    into  non- 
dimensional   form  by  using 

■f  =    -V  4>  W  =  We  u)  ,      r-ar'  (20) 


> 


where  W0  has  the  dimensions  of  a  velocity.   The  equa- 
tions are  then 


bb    5 

^ 


>r'    *t*  *V  J  V>    ^  ^  Ao   V        **'  I?) 


=  -r-'V,4^ 


(21) 


and 


J-(^i  >±>  _  <*±_  ^.\  -   C  -t    V.2^  <22) 


10 
where 

£    m    1*£a    ,    C  -  ££  (23) 


Now  ^   denotes  the  non-dimensional  form  of  the  operator. 
If  the  motion  is  slow  W0  may  be  considered  to  be  the 
velocity  component  along  the  center  of  the  pipe;  thus 
is  approximately  the  Reynolds  number.   The 


"V 

boundary  conditions  to  be  satisfied  are 


at  r'  *  1  .   Thus  for  laminar  flow  dynamical  similarity 
depends  solely  on  the  parameter  Kq 

Using  Goldstein'.^  3)  notation, 


This  may  be  replaced,  to  a  first  approximation,  by 


11 


where  Wm  is  the  mean  velocity  of  the  flow  in  the  curved 
pipe.   K  is  a  function  of  K  alone  since  W]n/W0  is  a 
function  of  K*  alone,  so  that  dynamical  similarity  de- 
pends on  K  .   (Dean  also  gave  the  equations  of  the  velocity 
profiles,  and  a  model  which  shows  the  profiles  has  been 
constructed.   See  Figure  2  .   The  secondary  flow  resulting 
from  the  curvature  is  indicated  by  the  streamlines  In  the 
plane  of  the  cross  section,  shown  in  Fig.  3  . ) 

Investigations  prior  to  Dean's  work  suffered  from 
lack  of  a  similarity  parameter  so  that,  though  it  was  known 
that  curvature  of  a  pipe  increased  the  resistance  to  flow 
above  that  of  a  straight  pipe,  the  results  of  experiments 
were  not  as  systematic  as  desired. 

Experiments  that  are  illustrative  of  these  early 
investigations  were:   those  of  Grindley  and  Gibson'^), 
conducted  to  determine  the  viscosity  of  air;  and  those  of 
Eustice^-5),  conducted  to  study  the  effects  of  curvature. 
The  tests  of  Grindley  and  Gibson  extended  over  the  range  of 
Reynolds  numbers  from  25  to  1^00  and,  hence,  did  not  cover 
the  transition  to  turbulence;  also  only  one  radius  of  cur- 
vature was  studied.   The  tests  of  Eustice  covered  a  wide 
range  of  Reynolds  numbers  (from  21  to  6000  for  one  pipe) 
and  curvatures,  but  the  pipes  did  not  have  a  circular 
cross  section.   Also,  no  Inlet  length  was  provided  for  full 
flow  development  ahead  of  the  testing  length. 

It  had  long  been  known  that  for  flow  In  a  straight 
pipe  the  resistance  is  proportional  to  the  velocity  of  flow 
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16 
when  the  flow  is  laminar,  but  it  ceases  to  be  proportional 
at  some  velocity,  which  is  known  to  be  the  critical  velocity 
at  which  transition  occurs.   It  was  also  known  that  curving 
a  pipe  caused  the  resistance  to  deviate  from  this  propor- 
tionality at  a  lower  velocity  than  that  of  the  straight 
pipe.   Early  investigators  thought  this  indicated  that 
turbulence  was  established  when  this  deviation  occurred,  so 
that  flow  in  a  curved  pipe  was  considered  to  be  less  stable 
than  that  in  a  straight  pipe. 

Dean,  however,  stated  that  the  secondary  flow  of 
the  curved  pipe,  which  is  the  cause  of  the  increase  in 
resistance,  is  a  steady  motion.   Eustlce  seemed  to  think  that 
there  apparently  was  no  critical  velocity  for  flow  in 
curved  pipes.   For  the  above  lack  of  better  experimental  data, 
Dean  seemed  to  accept  Eustlce' s  view. 

Shortly  after  the  publication  of  Dean's  paper, 
White^  '   conducted  resistance  measurements  on  pipes  of  three 
curvatures,  and  then  he  plotted  his  results  together  with 
those  of  Grindley  and  Gibson  as  a  function  of  Dean's  simi- 
larity parameter  and  apparently  confirmed  Dean's  analysis. 
Furthermore,  he  correctly  deduced  from  his  graphical  results 
that  transition  occurred  when  the  curve  for  a  particular 
curvature  ratio  left  the  main  curve,  i.e.,  ceased  to  be 
solely  a  function  of  Dean's  parameter.   The  graph  is  repro- 
duced as  Figure  4  .   For  a  curvature  ratio  of  1/2050 
transition  occurred  at  a  Reynolds  number  of  about  2300, 
while  a  ratio  of  1/15.5  required  a  Reynolds  number  of  about 
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9000.  Thus  it  appeared  that  the  minimum  critical  Reynolds 
number  for  flow  in  a  curved  pipe  depended  on  the  curvature 
ratio. 

That  this  was  so  was  at  once  verified  directly  by 
Taylor(?)  who  injected  a  stream  of  dye  into  the  flow  in  a 
helical  pipe  and  observed  the  onset  of  turbulence  visually. 
At  the  same  time  he  confirmed  Dean's  analysis  of  the 
secondary  flow  which  was  traced  out  by  Taylor's  color  band 
shown  graphically  in  Figure  5.   Taylor  writes, 


As  in  Reynolds'  experiments  with  a 
straight  pipe,  it  was  found  that  the  flow 
was  steady  up  to  a  certain  speed.  At  this 
speed  the  colour  band  began  to  vibrate  in 
an  irregular  manner,  but  it  still  seemed 
to  retain  its  identity  through  at  least 
one  whole  turn  of  the  helix.   This  indicates 
that  the  unsteadiness  was  not  at  first  of 
a  type  which  gives  rise  to  diffusion  of 
momentum  by  eddies  and  hence  to  a  rapid 
rise  in  resistance.   In  fig.  3,  which 
shows  graphically  the  results  of  these 
measurements,  the  point  where  this  irreg- 
ular vibration  first  made  Its  appearance 
is  Indicated  by  a  square,  thus   |5)  . 

On  increasing  the  speed  this  irreg- 
ular vibration  increased  in  violence  till, 
at  a  point  which  could  be  measured  with 
some  accuracy,  the  stream  became  completely 
turbulent  so  that  the  identity  of  the 
colour  band  was  lost  in  a  distance  of  a 
few  millimetres  from  Its  source.   The 
speeds  at  which  the  flows  became  completely 
turbulent  in  this  sense  are  shown  as 
circles,   ©  ,  in  fig.  3. 


(Taylor's  figure  3  is  reproduced  as  Figure  6.) 


Fig.  5«   Schematic  representation  of  Taylor's 

color  band  which  described  the  general 
type  of  motion  found  analytically  by 
Dean  for  laminar  flow. 
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Further , 

Inspection  of  fig.  3  shows  that  Kr. 
White's  conclusion  is  amply  verified  in 
experiments  with  colour  bands.  The  flow 
was  completely  turbulent  on  entering  the 
helix.  This  was  demonstrated  by  passing 
the  water  first  through  a  straight  glass 
tube  containing  a  1/2-mm.  hole  of  the 
same  type  as  that  in  the  glass  helix. 

At  approximately  the  speed  indicated 
by  Reynolds'  criterion,  ^P/^.  =   ^02-o  ,  the 
flow  became  turbulent  in  the  straight  pipe. 
It  then  entered  the  helix  in  this  state 
but  after  passing  through  only  one  coil  in 
the  case  of  glass  helix  No.  1,  or  2  1/2 
coils  in  the  case  of  helix  No.  2,  it  had 
become  steady  again.   This  steadiness 
persisted  up  to  a  value  dv/o^u.=  58  30   in 
the  case  of  helix  No.  1,  i.e.,  up  to  2.8 
times  the  speed  at  which  turbulence  persists 
in  a  straight  pipe.   It  will  be  noticed  in 
fig.  3  that  Kr.  White's  measurements  cor- 
respond far  more  closely  with  the  speed  at 
which  the  flow  first  appears  completely 
turbulent  than  with  the  speed  at  which  it 
first  begins  to  fluctuate.   This  might  be 
expected  because  a  fluctuation  which  does 
not  cause  the  colour  band  to  become  diffused 
is  not  likely  to  cause  any  considerable 
increase  in  resistance. 

Thus  White's  reasoning  was  justified  in  a  direct  manner. 


3*   The  Electric  Double  Layer 

Anderson  gave  a  rather  complete  account  of  the 
underlying  principles  of  the  electric  double  layer;  a 
brief  outline  of  this  is  given  here. 

From  studies  of  colloids  it  is  known  that,  at 
a  solid-liquid  interface,  the  solid  surface  becomes  charged. 
These  charges  are  supplied  by  the  liquid.   Since  there  is 
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25 
a  charge  on  the  interface,  there  must  also  be  an  equal 

and  opposite  layer  of  charge  adjacent  to  the  interface. 

These  charged  layers  are  called  the  electric  double  layer. 

Gouy^8)  and  Chapman^)  considered  the  double  layer  to  be 

ions  attached  to  the  wall  and  an  equivalent  ionic  charge 

of  opposite  sign  distributed  in  the  solution  near  the 

interface. 

The  charge  on  the  wall  is  treated  as  a  uniform 

surface  charge.   The  space  charge  in  the  solution  is  due 

to  the  unequal  distribution  of  ions,  assumed  to  have  no 

physical  size.   The  solvent  is  considered  as  a  continuous 

medium  entering  the  problem  only  through  its  dielectric 

constant. 

Assuming  an  Infinite  flat  wall,  the  differential 
equation  of  the  potential  is 


rt 


A% 


i=~^Zfen!0exP(-Z,ef/KT) 


(i) 


where  5      is  the  potential  (having  a  value  of    at 
the  wall ) , 

n±  is  the  number  density  of  the  1th  ion  species 

(having  a  value  of  nlo  at  the  wall), 

2T,  is  the  valence, 

€.  is  the  dielectric  constant, 

e  is  the  electronic  charge, 

k  Is  Boltzmann's  constant, 

and  T  Is  the  temperature  of  the  solution. 
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The  boundary  conditions  are 


d5 
S  »  O  ,   —  =  O         for    -X-^oo         (2) 

a  A 


where   X   is  the  distance  from  the  wall. 

-,*   4TT  e.  <_n.o  Z; 
Letting  X     ~    ■ (3) 


and  assuming  Z.e  yKT  to  be  small,  the  solution  ii 


or 


*  -  <  ^ 

5  -  50  e.  (4) 


showing  that  f   drops  by  a  factor  of  e"1   in  the  distance 

/x   •   This  value  is  called  the  thickness  of  the  double 
layer  or  Debye  length. 

**••   Streaming  Potential 

The  flow  of  the  fluid.  I.e.,  the  solvent,  sweeps 
out  ions  in  the  movable  part  of  the  double  layer,  developing 
a  streaming  potential  between  the  two  ends  of  the  passage. 
This  potential  difference  causes  a  conduction  current  to 
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flow  which  at  equilibrium  is  equal  and  opposite  to  the 
convection  current. 

For  laminar  flow  in  a  constant  area  pipe  caused 
by  a  pressure  gradient  which  is  constant  over  the  cross 
section,  analysis  yields 

AL  ^  i_L.  (5) 

where   E  is  the  potential  caused  by  the  fluid  flow, 

^u  is  the  absolute  viscosity, 

^  is  the  specific  conductance, 

and  p  is  the  pressure. 

Equation  (1)  shows  that  the  streaming  potential,  for  the 
assumptions  made,  is  a  linear  function  of  the  pressure 
drop  along  the  pipe. 

The  convection  current  is  given  approximately  by 


--M*.fcJ« 


where  vn     is  the  hydrodynamic  velocity  along  the  axis 
of  the  pipe, 

and  s      is  the  arc  length. 

Equation  (2)  shows  that  the  convection  current,  and  hence 
the  streaming  potential,  should  be  sensitive  to  the  velocity 
gradient  at  the  wall.   Thus  a  change  in  the  streaming 
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potential  should  provide  a  means  of  detecting  a  change 
dv» 


in 

dx 


5.   Results  of  Anderson's  Investigation 

Anderson  investigated  the  transition  to  turbulence 
in  a  straight  pipe  of  circular  cross  section  using 
fluctuations  in  the  streaming  potential  to  indicate  the 
onset  of  turbulence.   As  in  Taylor's  experiment,  it  was 
found  that  there  was  a  Reynolds  number  at  which  the  flow 
first  showed  some  instability  and  another  higher  Reynolds 
number  at  which  the  flow  was  completely  turbulent. 


II.   THE  EXFERIMENT 

1.  General  Description 

The  technique  developed  by  Anderson  for 
detecting  the  onset  of  turbulence  utilizes  the  very 
sharp  fluctuations  In  the  streaming  potential  which  are 
characteristic  of  turbulent  flow,  the  turbulent  eddies 
causing  changes  in  velocity  near  the  wall. 

2.  Laboratory 

The  laboratory  designed  by  Anderson  specifi- 
cally for  the  streaming  potential  measurements  obtains 
very  close  control  of  temperature  and  humidity.   The 
temperature  variations  must  be  low  because  of  the  strong 
dependence  of  the  fluid  viscosity  and  density  on  the 

temperature;  the  low  humidity  is  necessary  because  of  the 

-8      -Q 
very  small  electrical  currents  Involved:   10   to  10  7 

amps  with  a  source  leakage  on  the  order  of  50  megohms. 

The  temperature  of  the  room  can  be  controlled  to  ±0.75°F 

over  the  range  ^0  to  70°P,  which  is  a  continuous  variation 

over  a  b   minute  period.   The  humidity  of  the  room  can  be 

held  below  20%   R.H.   In  addition,  the  flow  system  was 

enclosed  in  a  Faraday  cage  to  avoid  unwanted  electrical 

noise. 
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3.   Apparatus 

A  schematic  drawing  of  the  flow  apparatus  Is 
shown  in  Pig.  7  .   The  reservoirs  R,  and  R2  are  5  liter 
Fyrex  boiling  flasks.   Tubes  L^  and  L2  are  6mm  Pyrex 
tubing.   A.,  and  A2  are  8mm  tubes  terminating  in  the 
sockets  of  ground  glass  ball  joints.   The  tube  T  is  the 
curved  pipe  which  is  a  precision  bore  capillary  tube. 
The  diameter  of  T  Is  ±  0.0002  inches.   The  length  of  T 
is  2k   inches.   The  ends  of  T  are  flared  into  8mm  tubes 
which  terminate  in  the  ball  portion  of  the  ball  joints. 
The  fluid  is  run  from  R.  to  R?  by  applying  gas  pressure 
through  V..   A  glass  trip  is  provided  at  B  to  ensure 
turbulent  flow  at  the  entrance  to  T.   Either  S1  or  S2  is 
connected  to  a  height  gage  to  determine  the  level  of  the 
fluid  in  the  flasks  to  within  t   0.002  inch. 

E-,  and  E~  are  Ag  -  AgCl  electrodes  based  on  a 
spiral  of  2k   AWG  platinum. 

Commercially  bottled  nitrogen  supplied  the  gas 
pressure.   A  schematic  diagram  of  the  pressure  system  is 
shown  in  Pig.   8  •   A  microscope  mounted  on  a  height  gage 
was  used  to  read  the  simple  mercury  manometer  (see  Pig.  9  ). 
With  this  system  and  the  accurate  bleed  valve  the  pressure 
could  be  controlled  to  within  t   0.005  inch  Hg. 


Fig.  9.   Manometer  and  optical  system  used  for 
the  pressure  measurements. 
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4.   Curving  and  Measuring:  of  the  Glass  Pipes 

The  curved  pipe  was  formed  from  a  precision 
bore  Pyrex  capillary  tube  which  was  curved  over  a 
precisely  machined  space  form  in  a  large  furnace. 
Approximately  10  minutes  in  the  furnace  at  a  temperature 
of  7^0°C  was  required  to  fabricate  the  capillary  tube. 
Any  greater  time  yields  an  unacceptable  distortion  of 
the  bore. 

The  diameter  of  the  pipe  in  the  plane  of  bending 
was  measured  at  3-lnch  Intervals  before  and  after  bending 
using  an  optical  comparator.   The  pipe  rested  in  a  tank 
of  glycerin  which  possesses  an  index  of  refraction  very 
close  to  that  of  Pyrex.   By  placing  the  bath  In  the 
temperature-controlled  laboratory  and  adjusting  the  temp- 
erature, the  index  of  refraction  of  the  glycerin  could  be 
brought  into  very  close  coincidence  with  that  of  the  Pyrex. 
This  occurred  when  the  outer  surface  of  the  glass  was  no 
longer  visible.   Tinted  glycerin  pumped  into  the  tube  made 
the  inner  surface  visible.   This  system  did  away  with  the 
troublesome  effects  of  refraction  from  the  glass  surfaces 
and  permitted  a  precise  measurement  of  the  tube  diameter. 

5«   Experimental  Procedures 

The  procedure  used  to  determine  the  critical 
Reynolds  number  was  to  measure  the  pressure  where  the 
turbulent  fluctuations  in  the  streaming  potential  first 
appeared.   This  pressure  was  then  used  to  find  the 


corresponding  Reynolds  number  from  a  calibration  curve 
obtained  from  mass  flow  measurements.  The  calibration 
curve  is  shown  in  Pig.  10  where 

P  is  the  pressure, 
d  is  the  diameter  of  the  tube, 
and  1   is  the  length  of  the  tube. 


CHAPTER  III 
RESULTS 

1.   1'easurements 

Using  the  procedure  previously  described  the 
critical  Reynolds  number  was  found  to  be  3190  with  a 
standard  deviation  of  20.   The  critical  point  was 
found  to  be  extremely  sensitive  to  temperature,  a  varia- 
tion of  1°C  resulting  in  a  change  in  the  critical 
pressure  of  about  1  inch  of  Hg. 

The  result  is  plotted  in  Fig.  11  as  the  point 
denoted  by  A  .      This  is  slightly  higher  than  the  values 
given  by  White,  which  were  obtained  from  mass  flow 
measurements.   This  apparent  difference  between  the  two 
techniques  is  attributed  to  the  need  for  longer  periods 
of  time  necessary  to  allow  for  sufficient  sampling  of 
the  fluctuations  in  the  streaming  potential  precisely 
at  the  critical  point.   Such  sampling  periods  are  not 
possible  with  the  present  system,  being  limited  by  the 
size  of  the  flasks. 
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APPENDIX   I 


The   Mavier-Stokes   Equation  for  a  Compressible 
Fluid   in  General   Orthogonal  Curvilinear  Coordinates. 
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